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ORIGINAL ARTICLE 

Epidemiological time series studies of PM 2 .5 
and daily mortality and hospital admissions: 
a systematic review and meta-analysis 

R W Atkinson, 1 S Kang, 1 H R Anderson, 1 ' 2 I C Mills, 3 H A Walton 2 ' 4 



ABSTRACT 

Background Short-term exposure to outdoor fine 
particulate matter (particles with a median aerodynamic 
diameter <2.5 jam (PM2.5)) air pollution has been 
associated with adverse health effects. Existing literature 
reviews have been limited in size and scope. 
Methods We conducted a comprehensive, systematic 
review and meta-analysis of 1 10 peer-reviewed time 
series studies indexed in medical databases to May 201 1 
to assess the evidence for associations between PM2.5 
and daily mortality and hospital admissions for a range 
of diseases and ages. We stratified our analyses by 
geographical region to determine the consistency of the 
evidence worldwide and investigated small study bias. 
Results Based upon 23 estimates for all-cause 
mortality, a 10 fxg/m 3 increment in PM 2 .5 was 
associated with a 1.04% (95% CI 0.52% to 1.56%) 
increase in the risk of death. Worldwide, there was 
substantial regional variation (0.25% to 2.08%). 
Associations for respiratory causes of death were larger 
than for cardiovascular causes, 1 .51 % (1 .01 % to 
2.01%) vs 0.84% (0.41% to 1.28%). Positive 
associations with mortality for most other causes of 
death and for cardiovascular and respiratory hospital 
admissions were also observed. We found evidence for 
small study bias in single-city mortality studies and in 
multicity studies of cardiovascular disease. 
Conclusions The consistency of the evidence for 
adverse health effects of short-term exposure to PM 2 . 5 
across a range of important health outcomes and 
diseases supports policy measures to control PM2.5 
concentrations. However, reasons for heterogeneity in 
effect estimates in different regions of the world require 
further investigation. Small study bias should also be 
considered in assessing and quantifying health risks 
from PM 2 . 5 . 



INTRODUCTION 

The adverse health effects of exposure to outdoor 
particulate matter air pollution are of concern to 
governments and health organisations worldwide. 1 2 
The evidence for these health effects has come 
from studies of the clinical, mechanistic and epi- 
demiological evidence of short-term and long-term 
exposures. While the epidemiological evidence 
relating short-term exposure to PMi 0 (particles 
with a median aerodynamic diameter <10 um) and 
related metrics (black smoke, total suspended parti- 
cles) with health effects is substantial, there are 
relatively few studies of fine particles measured as 
PM2.5 (particles with a median aerodynamic 



Key messages 



What is the key question? 

► Is there convincing and consistent evidence 
worldwide that short-term exposure to outdoor 
fine particulate matter (particles with a median 
aerodynamic diameter <2.5 urn (PM 2 . 5 )) air 
pollution is associated with increased risk of 
death and emergency admission to hospital? 

What is the bottom line? 

► We found evidence for adverse health effects of 
short-term exposure to PM 2 . 5 across a range of 
important health outcomes, diseases and age 
groups with substantial variation between 
different regions of the world that needs 
explanation. 

Why read on? 

► Our study provides a systematic, quantitative 
summary of the time series literature and 
reports new findings that suggest larger 
associations for respiratory causes of death 
than for cardiovascular causes and that the 
presence of publication bias in the literature 
could have important implications for public 
health policy. 

diameter <2.5 um). Reviews of the evidence 
linking exposure to PM2.5 to adverse health effects 
have relied upon a small number of published 
studies, restricted in health outcomes and geo- 
graphical coverage, or focused on differential 
PM2.5 toxicity. 3-11 

To summarise the existing evidence we conducted 
a comprehensive and systematic meta-analysis of 
time series studies of daily PM2.5 and daily mortality 
and hospital admissions published worldwide in the 
peer reviewed literature to May 2011. This included 
all disease outcomes for which there were sufficient 
studies for meta-analysis and combined results from 
single-city and multicity studies. We focused our 
analysis on single-pollutant rather than multipollu- 
tant models and upon all-year results in order to 
maximise the number of estimates available for 
inclusion in the review. We present estimates for 
WHO regions and assess between-region heterogen- 
eity. We also investigate whether there is evidence of 
publication (small study) bias between single-city 
study estimates and between multicity summary 
estimates. 
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METHODS 

Systematic ascertainment of relevant studies 

Time series (including case crossover) studies published in peer 
reviewed journals and indexed in online databases to May 2011 
(no start date specified) were identified via search strings using 
terms relating to study design, pollutant and health outcomes. A 
sifting process identified (from study titles, abstracts and the full 
paper) those time series studies suitable for inclusion in the 
review Study eligibility depended upon the details of the study 
design, statistical methods used and presentation of regression 
estimates and other data in numerical format. Further details are 
given in the online supplementary material. 

Extraction and coding of data 

Study details were entered into a Microsoft Access database 
(Microsoft Office 2010, Microsoft Corporation) and included 
citation information (title, authors, date of publication, etc) and 
details of effect estimates including health outcome (mortality 
or admission), diagnosis (International Classification of Diseases 
codes), age, and so on, and details of the pollutants (unit of 
measurement, range of exposure, etc). These data were used to 
calculate standardised effect estimates expressed as the percent- 
age change (and 95% CI) in the mean number of daily events 
associated with a 10 |xg/m 3 increase in PM 2 .5 concentration. 
The short-term relationships between air pollution and health 
effects are determined for given time lags (in days) between 
exposure and health events and investigators vary in which lags 
they study and report. 12 Hence an a priori lag selection protocol 
was devised and used to choose lag estimates for inclusion in 
the review without introducing bias (see online supplementary 
material for details). Additional data entry included the coding 
of the WHO region in which the study occurred (see online 
supplementary material, table SI). Studies were reviewed by a 
single statistician/epidemiologist before coding. All papers were 
read by RWA and data range checked prior to meta-analysis. 

Meta-analysis 

In the time series literature many cities have been studied on 
more than one occasion, in single-city studies and as part of a 
larger, coordinated multicity investigation. To ensure results 
from a city only appeared once in any one meta-analysis we 
applied an a priori estimate-selection protocol (see online 
supplementary material for details). 

We conducted meta-analyses only if there were 4+ single city 
estimates or if the set of estimates contained a multicity study 
summary estimate. Within each WHO region we conducted a 
two stage meta-analysis using a random-effects model for each 
stage. 13 In the first stage, single-city estimates were pooled to 
provide a summary estimate of the evidence from single-city 
studies. In the second stage, these summary estimates were 
pooled with the selected multicity study estimates to obtain a 
WHO region-specific summary estimate of the evidence. To 
assess heterogeneity between WHO regions we used the I 2 stat- 
istic which indicates the proportion of total variability between 
effect estimates due to heterogeneity. 14 I 2 statistics in the range 
0 to 30, 30 to 50 and >50 generally indicate low, moderate and 
high heterogeneity, respectively. Finally, a global summary esti- 
mate was calculated from WHO region-specific single-city 
summary estimates and multicity study estimates. 

Assessment of small study bias 

We investigated our selected single-city estimates and our 
pooled single-city and selected multicity estimates for evidence 



of small study bias using the methods of Begg and Egger. 15 16 
The former uses an adjusted rank correlation method to 
examine the association between study estimates and their vari- 
ance whereas the latter uses a regression approach. The impact 
of adjustment for small study bias was assessed using the 'trim 
and fill' method. 17 This method removes studies until symmetry 
in the funnel plot is achieved, recalculating the centre of the 
funnel before the removed studies are replaced together with 
their 'missing' mirror-image counterparts. A revised summary 
estimate is then calculated using all of the original studies, 
together with the hypothetical 'filled' studies. Our overall assess- 
ment of the evidence for small study bias was based upon the 
combined evidence presented by all three techniques. 

All analyses were conducted in STATA (STATA/SE V10, 
StataCorp Texas). 

RESULTS 

One hundred and ten time series studies of daily mortality (68) 
and hospital admissions (54) indexed in medical databases to 
May 2011, and providing numerical effect estimates, reported 
results for PM 2 .5 (see online supplementary material, table S2). 
Table 1 details the number of studies tabulated by outcome, 
disease, WHO region, age group and multicity versus single-city 
study design. The majority of studies of PM 2 .s and daily mortal- 
ity and hospital admissions have been conducted in North 
America and Europe with a small number of studies in other 
regions of the world. The most frequently reported estimates 
for daily mortality were for the all-ages group, followed by the 
65 + years group. For most populations the latter comprised a 
large proportion of the all-ages group so we confined our mor- 
tality analyses to the all-ages group. For hospital admissions we 
focused upon age-specific estimates in children and the elderly 
(ages 0-14 years and 65+ years, respectively). While the major- 
ity of studies were conducted in single cities, a substantial 
number reported findings from multicity studies. 

Mortality 

Summary estimates (95% CIs) per 10 |Jig/m 3 increment in PM 2 .5 
and all-age, all-cause and cause-specific mortality are presented 
in figure 1. All associations were positive and for all, except 
chronic COPD, lower CIs were above unity. For all-cause mor- 
tality, 23 single-city and multicity study estimates were selected 
for meta-analysis from the 43 estimates identified in the review 
(see online supplementary material, figure SI). The overall 
random effects summary estimate was 1.04% (95% CI 0.52% 
to 1.56%) per 10 |xg/m 3 increment in PM2.5. WHO region- 
specific summary estimates varied substantially (I 2 =93%) from 
0.25% to 2.08% (table 2). 

While fewer estimates were available for cardiovascular (see 
online supplementary material, figure S2) and respiratory (see 
online supplementary material, figure S3) mortality, the overall 
summary estimate for all respiratory causes of death was larger 
than for all cardiovascular causes, 1.51% (95% CI 1.01% to 
2.01%) versus 0.84% (95% CI 0.41% to 1.28%), respectively. 
For both causes of death, associations were positive in all WHO 
regions (table 2) and heterogeneous for cardiovascular deaths 
(I 2 = 76%) but not respiratory deaths (I 2 = 0%). Associations 
between PM2.5 and death from ischaemic heart disease, stroke 
and COPD were 3.36% (0.68%, 6.10%), 1.85% (0.74%, 
2.97%) and 2.86% (-0.12%, 5.93%) per 10 |JLg/m 3 , respect- 
ively, although the evidence was restricted to a small number of 
single-city and multicity studies (see online supplementary 
material, table S3 and figures S4-S6). 
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Table 1 Time series studies of PM 2 . 5 and mortality and hospital admissions 



Multicity study 10131 Multicity study Single-city study 



Outcome 




Mortality 


Hospital admission 


Mortality 


Hospital admission 


Mortality 


Hospital admission 


Total 




68 


54 


17 


11 


51 


43 


Disease 


Respiratory 


33 


43 


7 


9 


26 


34 




Cardiovascular 


41 


34 


9 


9 


32 


25 




All-cause 


56 


2 


15 


0 


41 


2 




Other 


7 


3 


2 


2 


5 


1 


WHO region 


American region A 


33 


31 


13 


8 


20 


23 




European region A 


20 


10 


2 


1 


18 


9 




Western Pacific region B 


6 


6 


0 


0 


6 


6 




American region B 


6 


2 


0 


0 


6 


2 




Western Pacific region A 


4 


4 


3 


2 


1 


2 




South-East Asia region D 


0 


1 


0 


0 


0 


1 


Age group 


All ages 


54 


21 


16 


1 


40 


20 




Elderly 


26 


28 


5 


9 


21 


19 




Not elderly 


4 


4 


1 


1 


3 


3 




Adult 


1 


2 


0 


0 


1 


2 




Young adult 


0 


9 


0 


2 


0 


7 




Children 


1 


18 


0 


3 


1 


15 




Other 


2 


3 


0 


0 


2 


3 



Hospital admissions 

Table 3 gives summary estimates for all-age, cardiovascular and 
respiratory causes of hospital admissions with individual study 
results presented in online supplementary material, figures S7 
and S14. 

PM 2 .5 concentrations were positively associated with increases 
in risk of admission for cardiovascular diseases, 0.90% (95% CI 
0.26% to 1.53%) and respiratory diseases, 0.96% (95% CI 



-0.63% to 2.58%) per 10 |JLg/m 3 , respectively, with heterogen- 
eity between WHO regions for respiratory diseases (I 2 = 80%) 
but not cardiovascular diseases (I 2 =0%). Figure 2 illustrates 
summary estimates for specific cardiovascular diseases in ages 
65+ years and for respiratory diseases in ages 65+ years and 
children aged 0-14 years. All associations were positive except 
for stroke and for all, except COPD including asthma, lower CIs 
exceeded 0%. Details of WHO-speciflc summary estimates are 



All cause - 



Cardiovascular- 



IHD- 



Stroke 



Respiratory - 



COPD 



i- 



i — •- 



1 1 1 1 

-2 0 2 4 6 8 
% (95% CI) per 10 ug/m3 increase in PM2.5 



Figure 1 Summary estimates (95% confidence intervals) for all-cause 
and cause-specific mortality. 



Table 2 Meta-analysis results for all-age, all-cause mortality and 
cause-specific mortality by WHO region 



WHO region 


All* 

(SC/MC) 


Selectedt 
(SC/MC) 


RE (95% CI)* 


I 2 

(%)§ 


All Cause 










AMR A 


13/12 


5/2 


0.94 (0.73 to 1.16) 


93 


AMR B 


4/0 


2/0 


2.08 (1.60 to 2.56) 




EUR A 


12/1 


9/1 


1.23 (0.45 to 2.01) 




WPR A 


0/1 


0/1 


0.90 (-0.70 to 2.53) 




WPR B 


5/0 


3/0 


0.25 (0.06 to 0.44) 




Summary H 




4/4 


1.04 (0.52 to 1.56) 




Cardiovascular 










AMR A 


10/3 


6/1 


0.84 (0.47 to 1.20) 


76 


AMR B 


3/0 


2/0 


0.13 (-0.71 to 0.98) 




EUR A 


6/1 


6/1 


2.26 (1.23 to 3.29) 




WPR B 


4/0 


2/0 


0.56 (0.31 to 0.81) 




Summary H 




4/2 


0.84 (0.41 to 1.28) 




Respiratory 










AMR A 


4/5 


4/1 


1.39 (0.62 to 2.16) 


0 


AMR B 


3/0 


2/0 


0.88 (-1.88 to 3.71) 




EUR A 


7/0 


7/0 


3.81 (0.57 to 7.16) 




WPR B 


4/0 


2/0 


1.49 (0.04 to 2.96) 




SummaryU 




4/1 


1.51 (1.01 to 2.01) 





*Numbers of single-city(SC)/multicity (MC) estimates available from all studies. 

tNumbers of single-city(SC)/multicity (MC) estimates selected for meta-analysis (see 

estimate selection protocol in Methods section). 

^Random effects summary estimate (95% CI) per 10 (ig/m 3 . 

§l 2 statistic for heterogeneity. 

HEstimate numbers for 'Summary' refers to the number of pooled (from single-city 
estimates) and multicity estimates used to calculate the overall summary estimate 
across WHO regions. 

AMR, Region of the Americas; EUR, European Region; WPR/SEAR, South East Asian Region. 
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Table 3 Meta-analysis results for all-age, cardiovascular and 
respiratory hospital admissions by WHO region 





All* 


Selectedt 




I 2 


WHO region 


(SC/MC) 


(SC/MC 


RE (95% CI)* 


(%)§ 


Cardiovascular 










AMR A 


2/0 


1/0 


0.00 (-2.85 to 2.93) 


0 


EUR A 


4/1 


4/1 


0.91 (0.17 to 1.66) 




WPR A 


1/0 


1/0 


1.04 (-0.30 to 2.39) 




Summaryll 




3/1 


0.90 (0.26 to 1.53) 




Respiratory 










AMR A 


1/0 


1/0 


-2.00 (-6.00 to 2.17) 


80 


EUR A 


3/0 


3/0 


1.90 (-0.18 to 4.02) 




SEAR D 


1/0 


1/0 


0.12 (0.08 to 0.16) 




WPR A 


1/0 


1/0 


2.38 (1.04 to 3.73) 




Summaryll 




4/0 


0.96 (-0.63 to 2.58) 





*Numbers of single-city(SC)/multicity (MC) estimates available from all studies. 

tNumbers of single-city(SC)/multicity (MC) estimates selected for meta-analysis (see 

estimate selection protocol in Methods section). 

{Random effects summary estimate (95% CI) per 10 (ig/m 3 . 

§l 2 statistic for heterogeneity. 

HEstimate numbers for 'Summary' refers to the number of pooled (from single-city 
estimates) and multicity estimates used to calculate the overall summary estimate 
across WHO regions. 



given in online supplementary material, tables S4 and S5, and 
study-specific estimates in online supplementary material, figures 
S7-S13 and S15-S20. 



Small study bias 

The impact of adjustments for small study bias in single-city 
studies within WHO regions and between the pooled single-city 
and multicity estimates for all-age, all-cause and cause-specific 
mortality are shown in table 4. We found evidence for small 
study bias in single-city mortality studies and in multicity studies 
of cardiovascular disease. 

Figure 3 illustrates, for cardiovascular mortality, this adjust- 
ment using a funnel plot of individual study estimates and 
showing the random effects summary estimates with, and 
without, inclusion of the two 'filled' estimates identified by the 
'trim and fill' procedure. We did not find evidence of small 
study bias in either single-city or multicity estimates for cardio- 
vascular or respiratory hospital admissions (data not shown). 

DISCUSSION 

In this systematic review we identified 110 ecological time series 
studies of short-term exposure to outdoor PM 2 .5 and daily mor- 
tality and hospital admissions indexed in medical databases up 
to May 2011. Our meta-analysis of effect estimates from these 
studies indicated positive associations with daily all-cause and 
cause-specific mortality and cause-specific and age-specific hos- 
pital admissions, with some evidence of heterogeneity. We also 
found evidence for small study bias in single-city estimates and 
between pooled and multicity estimates. 

There are a number of plausible biomedical explanations for 
associations between short-term exposure to fine particles and 
adverse health outcomes. 5 6 It is hypothesised that small effects 
cause clinical events when experienced by individuals who are 
already vulnerable due to existing chronic or acute disease. Our 
review indicates that such effects are observed even at the rela- 
tively low levels of fine particles found in developed countries. 
Our results reinforce the public health importance of fine parti- 
cles on health. While the estimates are small, the impact is sub- 
stantial because the entire population is exposed. Impact 



Cardiovascular, All Ages 
Cardiovascular, 65+ yrs 
IHD, 65+ yrs 
Heart Failure, 65+ yrs 
Cardiac, 65+ yrs 
Stroke, 65+ yrs- 
Dysrhythmia, 65+ yrs 
Respiratory, All Ages - 
Respiratory, 65+ yrs- 
COPD incl asthma, 65+ yrs- 
COPD excl asthma, 65+ yrs- 
Lower Resp Infection, 65+ yrs- 
Respiratory, 0-14 yrs- 
Asthma, 0-14 yrs- 



i — • — i 
i •— 



I • 1 



-• 1 



I • 1 



— I 1 1 1 

-2 0 2 4 6 8 

% (95% CI) per 10 ug/m3 PM2.5 



Figure 2 Summary estimates (95% confidence intervals) for 
cardiovascular and respiratory hospital admissions. 



assessments of PM 2 .5 on mortality are based on cohort rather 
than time-series evidence because this enables years of life lost 
to be estimated. 18 19 However, most cohort evidence is from 
North America or Western Europe. Our finding that short-term 
associations occur worldwide supports the generalisation of 
cohort based estimates globally 19 while at the same time indicat- 
ing that there may also be some heterogeneity. 

Our study extends the literature reviewing the evidence for 
health effects of short-term exposure to PM 2 .5 derived from 
time series studies. The numbers of single-city and multicity 
studies have increased substantially in recent years, from 55 in 
2005 11 to 110 identified in this review. In 2009 the 
Environmental Protection Agency (EPA) summarised evidence 
from studies of mainly US populations 5 published between 
2002 and May 2009. Also in 2009, the American Heart 
Association issued an updated statement on the effects of par- 
ticulate air pollution on cardiovascular disease including evi- 
dence from time series studies published to March 2009. 6 A 
more recent review in 20 12 8 focused upon studies of PM 2 .5 
components indexed in the Science Citation Database up to 
October 2010. Our study complements these previous reviews 
by providing a meta-analysis for a much larger, more recent 
(indexed in medical databases to May 2011) and broader litera- 
ture incorporating studies irrespective of geographical location 
and disease outcome. 

Across the five WHO regions studied, our summary estimates 
for all-age, all-cause mortality ranged from 0.25% to 2.08% 
with an overall estimate of 1.04%, comparable with the 
summary estimate of 1.2% derived from seven studies by Levy 
et al. 8 While we found evidence of statistical heterogeneity 
across region-specific estimates both values are consistent with a 
mortality hazard from short-term exposure to fine particulate 
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Table 4 Assessment of bias in single-city studies and in pooled estimates for all-cause mortality and cause-specific mortality 



All-cause Cardiovascular Respiratory 



No Adjustment* 


1.04 (0.52 to 1.56) 


0.84 (0.41 to 1.28) 






1.51 (1.01 to 2.01) 


Single-city biast 












WHO region 


Amr A 


Amr A 


Amr B 


Eur A 


Amr B 


p Value Begg 


0.14 


0.57 


0.32 


0.35 


0.32 


p Value Egger 


0.003 


0.42 


NA 


0.32 


NA 


# Estimates 


5 


6 


2 


6 


2 


#Trim and fill 


8 


7 


3 


7 


3 


Single-city* 


0.97 (0.46 to 1.48) 


0.78 (0.35 to 1.21) 






1.00 (0.73 to 1.27) 


Multicity bias§ 












p Value Begg 


0.81 


0.09 






0.46 


p Value Egger 


0.36 


0.32 






0.52 


# Estimates 


8 


6 






8 


# Trim and fill 


8 


8 






8 


Pooled single and multicityll 


0.97 (0.46 to 1.48) 


0.57 (0.09 to 1.05) 






1.00 (0.73 to 1.27) 



*Random effects summary estimate (95% CI) per 10 jig/m 3 without adjustment for small study bias. 

tAnalysis of bias in single-city studies by WHO region (where found): Begg's test p value, Egger's test p value, number of estimates prior to application of 'trim and fill' technique, 
number of estimates after application of 'trim and fill' technique. 

^Overall summary estimate calculated after application of 'trim and fill' technique to single-city estimates by WHO region. 

§Bias between pooled single-city estimates and multicity estimates, Begg's test p value, Egger's test p value, number of estimates prior to application of the trim and fill' technique, 
number of estimates after application of the 'trim and fill' technique. 

HOverall summary estimate after application of the 'trim and fill' technique to single-city estimates within WHO region and between pooled single-city estimates and multicity estimates. 



matter. The reasons for this variability in effect estimates 
between WHO regions warrant further investigation but may 
reflect variations in population vulnerability and/or differential 
toxicity of sources, pollutant mixtures and pollution 
monitoring. 

Our finding that the association for respiratory mortality 
(1.51%) was larger than for cardiovascular mortality (0.84%) — 
a finding observed within all WHO regions — was also reported 
by the EPA in their review. 5 However, the EPA noted the coher- 
ence in associations between PM 2 .s and cardiovascular mortality 
and morbidity outcomes and the lack of such coherence for 
respiratory diseases. We were only able to assess the extent of 
coherence by comparing associations for mortality and admis- 
sions. For cardiovascular diseases our overall summary estimates 
for PM2.5 and admissions and mortality were comparable con- 
firming the coherence reported by the EPA although we note 
substantial disagreement between region-specific summary 



♦I 



-10 12 3 

(3x1000 

Figure 3 Funnel plot of pooled single-city and multi-city summary 
estimates for cardiovascular mortality including 'filled' estimates. 
Random effects summary estimates without (long-dash line) and with 
(short-dash line) adjustment using the Trim & Fill procedure. 



estimates for each outcome including North America. We also 
note the negative, although statistically not significant, associa- 
tions observed in three single-city studies in Europe and the 
Western Pacific region for hospital admissions for stroke in 
adults over the age of 65 years. For respiratory diseases, we 
found that the overall summary estimates for PM2.5 and mortal- 
ity and admissions were broadly comparable although for 
respiratory admissions the summary estimate was approximately 
two-thirds that for mortality and the CI for the summary esti- 
mate straddled 0% due to the negative association reported in a 
single-city study from North America. We note however that 
while the same ICD codes are used for mortality and admis- 
sions, the way they are used is different — underlying cause of 
death for the former and immediate cause of admission for the 
latter. This might affect comparability of certain categories such 
as pneumonia or heart failure though they would still fall within 
the broad cardiovascular or respiratory rubric. 

The main strengths of our study are: (1) a protocol driven 
approach to the identification, coding and selection of effect 
estimates for meta-analysis to minimise selection bias through- 
out the review process; (2) inclusion of all health outcomes for 
which sufficient estimates were available for meta-analysis; (3) 
no limitations on study location or language, and (4) stratifica- 
tion of results by WHO region. However, in common with 
other reviews, our study is limited by: (1) the need for numer- 
ical, rather than graphical, presentation of data to facilitate 
quantitative meta-analysis; (2) the authors' choice of results to 
publish; and (3) having no assessment of the 'grey' literature. It 
is therefore important to assess the extent to which the infer- 
ences and quantitative estimates presented here have been 
affected by small study bias, a process that leads to the published 
literature being unrepresentative of the totality of evidence. 20 

Our analysis of small study effects, a generic term also encom- 
passing publication bias and a range of other potential biases, 21 
suggests that this may be an important issue in the meta-analysis 
of single-city and multicity estimates. The former has been 
noted previously in relation to time series studies of PMk) 11 but 
the observation that meta-analysis of multicity estimates can be 
similarly affected is new. Sterne et al suggested that the greater 
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investment of time and money in larger studies (such as multi- 
city time-series studies) meant such studies would be more likely 
to be of high methodological quality and published even if their 
results are negative. 21 Sterne et al also noted that the 'trim and 
fill' procedure detects 'missing' studies in a substantial propor- 
tion of meta-analyses, even in the absence of bias. Thus, there is 
a danger that the application of the procedure could mean 
adding and adjusting for non-existent studies in response to 
funnel plot asymmetry arising from nothing more than random 
variation. Our finding that the 'trim and fill' adjustment substan- 
tially reduced the magnitude and precision of the associations 
between PM2.5 and cardiovascular mortality should therefore be 
interpreted with some caution, especially so given the unre- 
markable p values from the Begg and Egger tests and the funnel 
plot presented in figure 3. 

Other potential sources of bias in our study methodology 
could be: (1) the selection of estimates from the individual 
papers; and (2) the selection of study estimates for meta-analysis. 
Bias arising from the former is possible given the tendency of 
investigators to assess associations at different time lags between 
exposure and health event. To address this we used a protocol for 
estimate selection that was independent of the direction of an 
association. In previous work this has been shown not to intro- 
duce bias but increase between-estimate heterogeneity 11 so we 
believe that it is unlikely that we are overstating the magnitude of 
associations. Our protocol for the selection of estimates (from 
those available) for meta-analysis similarly did not consider direc- 
tion of associations, instead basing the selection upon geograph- 
ical coverage, publication date and length of study period. It also 
ensured that no single location appeared more than once in a 
summary estimate. While our approach was just one of many ana- 
lytical strategies that could have been adopted, we believe that 
the resulting summary estimates were unlikely to be systematically 
biased since none of our selection criteria included the direction 
or magnitude, of the individual study effect estimates. The rele- 
vance of these issues cannot be overstated since results from 
meta-analyses can be important in underpinning the use of more 
limited cohort data in worldwide health impact assessments. 19 

Our review points to adverse associations between short-term 
exposure to daily concentrations of PM2.5 and daily mortality 
and hospital admissions across a range of diseases and age groups 
which supports continued policy measures to control PM2.5 
levels worldwide. However, we note that the evidence for these 
associations is concentrated in a small number of geographical 
regions of the world and also limited to the broader categories of 
disease. Further studies from other developed countries, in par- 
ticular Asia and Eastern Europe, are needed to confirm the 
observed associations. Also, new studies including specific, rather 
than broad, categories of diseases would increase understanding 
of the populations at risk and may also add to our understanding 
of mechanism of effect. In addition, the reasons for the hetero- 
geneity in effect estimates in different regions of the world 
require further research as they may be relevant to the formula- 
tion of policy measures. Our findings also suggest that, for the 
purpose of health impact assessment, some consideration of, and 
adjustment for, small study bias in results from multicity as well 
as single-city studies should be undertaken. 
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